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Abstract 

The decay of a parent particle into two or more daughter particles results in an entangled 
quantum state, as a consequence of conservation laws in the decay process. We use the Wigner- 
Weisskopf formalism to construct an approximation to this state that evolves in time in a manifestly 
Qh' unitary way. We then construct the entanglement entropy for one of the daughter particles by use 

■^ . of the reduced density matrix obtained by tracing out the unobserved states and follow its time 

CN ! evolution. We find that it grows over a time scale determined by the lifetime of the parent particle 

to a maximum, which when the width of the parent particle is narrow, describes the phase space 
distribution of maximally entangled Bell- like states. 
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I. INTRODUCTION 

Once described as the source of "spooky action at a distance" by Einstein, Podolsky and 
Rosen (EPR) |l|, quantum entanglement has now come to viewed as a resource to be exploited 
in a number of venues. It serves as the workhorse for quantum computations (2H8|, as well as 
a way to measure a variety of quantities in particle physics, such as flavor entanglement in 
T(4S') — )■ B'^B via the analysis of the time dependence of the flavor asymmetry[9l] and time 
reversal violation in the B^ system [lOl. Ill| . The possibility has also been advanced 12|, [l3 



to use entanglement correlations to establish bounds on the Bs — Bg width difference and 
CP violating phases. Entanglement between the charged lepton and its associated neutrino 
in the decay of pseudoscalar mesons has been recently argued to pla y a n important role 



in the coherence (and decoherence) aspects of neutrino oscillations |14j-|l 7l| with potentially 



important corrections in short baseline oscillation experiments |18l]. In condensed matter and 
quantum optics, the spontaneous decay of excited atomic states leads to quantum entangled 
states of photons and atoms or spin-qubits, which can then be exploited as platforms for 
quantum computingjsl, |j|. 

Given a pure quantum system consisting of entangled subsystems, it may not be possible 
to measure the separate state of all of the subsystems (or while possible, we may opt not 
to measure them). We can then construct a reduced density matrix for the subsystem(s) we 
do measure by tracing over the allowed states of the unobserved subsystem(s). This then 
leads directly to the concept of entanglement entropy: this is the von Neumann entropy of 
the reduced density matrix. It reflects the loss of information that was originally present 
in the entangled state from the quantum correlations. This entanglement entropy has been 
the focus of several studies in statistical and quantum field theories where subsystems are 
spatially correlated across boundaries by tracing over the degrees of freedom of one part of 
the system pj|-|2J] and have been extended to the case including black holesj25|, |26| , particle 



production in time dependent backgrounds [2 7l| and cosmological space tim es [28l| . Momentum 
space entanglement and renormalization has been recently studied in ref. [2i 

While entanglement entropy has been mostly studied within the context of a quantum 
system subdivided by space-like regions (see references above), in this article we study the 
time evolution of the entanglement entropy in the ubiquitous case of particle decay. The 
construction of the relevant states in this case relies on the Wigner-Weisskopf theory of 



spontaneous emission '3d. l3l[|. which provides a non-perturbative method for obtaining the 
quantum state arising from spontaneous decay. The knowledge of the full quantum state 
can then be used to obtain the entanglement entropy contained in, for example, the photon- 
spin qubit correlations generated from the dynamics of spontaneous decay in solid state 
systems J32l]. 

Recently this theory was generalized to relativistic quantum field theory to yield insight 
into the quantum states from particle decay in cosmology [33|, and to describe potential 
decoherence effects in neutrino oscillations in short baseline experiments [1 71. Il8| . 

Here we extend and generalize the Wigner-Weisskopf method discussed in|l8l. 1321 . |33| to 
describe particle decay in quantum field theory and apply it to the simple case of a bosonic 
parent particle decaying into two bosonic daughter particles, although we argue that the 
results are general. We address in detail the important aspect of unitarity and obtain the 
entanglement entropy by tracing over the degrees of freedom associated with an unobserved 
daughter particle. We show in detail how unitary time evolution yields an entanglement 
entropy that grows over the lifetime of the parent particle and saturates to the entropy of 



maximally entangled states. 

II. THE WIGNER-WEISSKOPF APPROXIMATION 

Consider a system described by a total Hamiltonian that can be decomposed as if = 
Hq + Hi, where Hq is the free Hamiltonian and H^ is the interaction part. As usual, the 
time evolution of a state in the interaction picture is given by 

^jtMt)) = Hj{t)Mt)) (ii.i) 

where 

Hi{t) = e'^'^'Hie-'^"' (11.2) 

is the interaction Hamiltonian in the interaction picture. The formal solution of (]11.1|) is 
given by 

|M/(t)) = f/(t,to)|^(to)) ; f/(t,to) = T(e~^^o^^(*')'^*') (H.3) 

Expanding the state |\I'(t)) in the basis of eigenstates of Hq we have 

mt)) = J2^n{t)\n) (11.4) 



where HqIti) = En\n). From eq.f lll.ll) . upon expanding in basis states \n), it follows that 



^^ = -tJ2{^\Hi{t)\m)Cm{t). (11.5) 

m 

This is an infinite set of differential equations that can be solved hierarchically leading 
in general to integro-differential equations at any given order in the perturbative expansion. 
Progress can be made by considering the evolution of the states that are coupled at a given 
order in perturbation theory and solving the coupled equations for these, thereby truncating 
the hierarchy at a given order in the perturbative expansion. 

For the case of interest here, consider an initial value problem in which the system is 
prepared at an initial time t = in a state \A), so that Ca{0) = 1, Cn^A = and that the 
Hamiltonian Hj couples these states to a set of states \k,). Then, we can close the hierarchy 
at second order in the interaction by keeping only the coupling of the states \A) ^ \k), i.e. 

^CA{t) = -^J2{^\HJ{t)\K)C^{t) (11.6) 



dt 
d 



^C^{t) = -i{K\Hi{t)\A)CA{t). (11.7) 

dt 



Using the initial conditions we obtain 



a(t) = -i f dt'{K\Hi{t')\A)CA{t') , (11.8) 

Jo 



which when inserted into (111.6^ leads to 

jCA{t) = - j dt'Y,{MHi{t)\K){K\Hi{t')\A)CA{t') = - j dt'J:Ait-t')CAit'), (II.9) 
where the second order self energy has been introduced 
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j:^(t-t')^J2^A\Hjit)\K){K\Hjit')\A) = J2\{A\HiiOM e^(^-"^«)(*-*') . (11.10) 

Higher order corrections can be included by enlarging the hierarchy, i.e. by considering 
the equations that couple the states k to other states k' via the Hamiltonian. The coefficients 
for the states n' can be obtained by integration and can be inserted back in the equations for 
the coefficients C^ (which already include their coupling to \A)). Then formally integrating 
the equation and inserting the results back into the equation for 1^4) generates higher order 
corrections to the self-energy S^. Finally, solving for the time evolution of Ca allows us to 
obtain the time evolution of the other coefficients. 

As seen by the procedure described above, the Wigner-Weisskopf approach can be used 
to construct an approximate version of the quantum state in the presence of interactions. 
However, what is not altogether obvious is that the truncation of states used to construct the 
state gives rise to a state whose time evolution is unitary. This will be extremely important 
in the sequel since we will want to follow that time evolution of the entanglement entropy 
this state would provide after tracing out an unobserved subsystem as discussed in the 
introduction. We will need to be sure that there are no spurious effects in this evolution due 
to an approximation to the state. 

The statement of unitary is one of conservation of probability. From the evolution equa- 
tion flll.Sp and its complex conjugate it follows that 



tX^cm?--'Y. 



dt 



m,n 



CUt)C:it){n\Hi{t)\m)-Cn{t)C*^{t){m\Hj{t)\n) 



0, (11.11) 



as can be seen by relabeling m ^ n in the second term. Therefore Yin l^™(^)P ~ constant. 
Now this is an exact result; the question is whether and how is it fulffiled in the Wigner- 
Weisskopf approximation obtained by truncating the hierarchy to the set of equations 

fiii.6yii.7p . 

Using (IIL81[ILT0|1 ) consider 

$^|a(t)|' = / dt^C\{t,) [ dt2^A{ti,t2)CA{t2) (11.12) 

inserting 1 = 6(ti — ^2) + ©(1^2 — ti) in the time integrals it follows that 

V|a(t)|2 = / rftiQ(ti) r dt2T.A{t,M)CA{t2) 

+ f dt2CA{t2) r dtj:A{tiMC\{t{) (11.13) 

JO JO 



using T,A{ti,t2) = S^(t2,ti) relabelling ti -H- ^2 in the second line of (111.13^ and using (1II.9P 
we find 



J2\^M' 



dti 



C\{U)^CA{h)+CA{U)^C\{U) 



d 



,Q ail 

(11.14) 
where we have used (Ca(0) = 1). This is the statement of unitary time evolution, namely 



ic^A(t)r+$^ia(t)r = ic^ 



(11.15) 



This is an important result. As we will see below, standard perturbative calculations 
of this state would not yield a state that evolved unitarily. The Wigner-Weisskopf state 
involves a non-perturbative dressing up of the state and despite its approximate nature, this 
dressing up captures the physics to a sufficient extent to guarantee unitary time evolution. 

III. REDUCED DENSITY MATRIX AND ENTANGLEMENT ENTROPY. 



Now we are ready to turn to the problem we really want to consider: the state that 
appears after the decay of a parent particle $ into two daughters x-, '4'- For simplicity and to 
highlight the main concepts we treat all fields as bosonic massive fields with masses m<j, and 
m^, m^ respectively. We consider a typical interaction vertex described by the interaction 
Hamiltonian 

Hi = g I d^x $(x) x{.x) ^(x) . (III.l) 



We quantize the fields in a volume V , so that in the interaction picture they can be 
written as: 



^(f,t) = ^ 



vmy 



a^g-iSfeigifc-5'^^^giE.tg-ifc.: 



^ = *,X,V^- 



(III.2) 



For the case of interest here, namely the decay process $ — )■ X'?/' we consider that the 
initial state is a single particle $ at rest, therefore the initial condition is (7$ (A; = 0; t = 0) = 
1 , Ck = for \n) 7^ I if)- The interaction Hamiltonian ( IIII.ip connects the initial state, 
|1?) to the states, \n) = \x-p) ® li^p)- These states in turn are coupled back to |1?) via Hj, 
these processes are depicted in fig. (II]). 



X 



$ 




X 



(^ 



FIG. 1: Transitions \^) O |x)l^) up to order g^ that determine S$. 



Thus to leading order in g we find from tlie intermediate states sliown in fig. ([I]) 



E^it-f) = j2\{^l\Himx~p,i^p)\ 



2 i(m^-E^{p)-E^(p))(t-t') 



(III.3) 



The interaction Hamiltonian also connects a single $-particle state to an intermediate 
state with three other particles and this state back to the single $ particle state yielding a 
disconnected contribution to the self energy depiced in fig.©. This contribution is just a 
renormalization of the vacuum energy and only contributes to an overall phase that multiplies 
the single particle $ state and will be neglected in the following analysis. For a more detailed 
discussion of this contribution see ref. 33 . 
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FIG. 2: Order g^ correction to the vacuum energy. Yields an overall phase for the quantum state 



At higher order in g there are higher order contributions to the $ self energy from other 
multiparticle states. However, we will only be considering states connected to |1*) to first 



order in perturbation theory. As discussed in detail in refs. [18|, |33| the Wigner-Weisskopf 
method provides a non-perturbative resummation of self-energies in real time akin to the 
dynamical renormalization group |34l| . The self-energy eq.f lIII.3|) is recognized as the one- 
loop retarded self energy of the field $ with the |x) 1"^) intermediate state [l8|, |33| . 

Solving eq. flll.Qp produces a solution for the time evolution of the $ amplitude. We can 
use the solution for C$(t) to obtain an expression for the amplitudes Ct^{t) which allows for 
computation of the probability of occupying a particular state at any given time. We may 
solve eq. ( (]II.9|) ) either via Laplace transform, or in the case of weak coupling, a derivative 
expansion which yields the same result at long times (t ^ l/niq,). In ref. 33| the equivalence 



between the two approaches is discussed in detail. Here, we follow the latter method which 



is the original Wigner-Weisskopf approximation J30|, |31 
We begin by defining the quantity 



Wo{t,t') 



so that 



dt"S$(t-t") 



(III.4) 



d 



— I^o(t,t') = S-,(t-t') , W^o(t,0)=0 



dt 



(IH.5) 



Integrating eq.( ]II.9p by parts yields 

jC^{t) = - I dt'J:^{t~t') C^{t') = -Wo{t,t) C$(t) + f dt'Wo{t,t')-^C^{t'). (III.6) 

The first is term second order in Hj whereas the second term is of fourth order in Hj 
and will be neglected. This approximation is equivalent to the Dyson resummation of the 
one- loop self energy diagrams. Thus to leading order, eq. ( 111.91) becomes 



j^C^{t) + Wo{t,t)C^{t) = 0, (III.7) 



where 



Wo{t,t)= [ dt'T.^{t - t') = [ rft'V|(l?|^,(0)|x-p,^p)|V('"*-^^(P)-^*(p))(*-*') (111.8) 
^0 Jo p. 

Inserting a convergence factor and taking the limit t — ?■ oo consistently with the Wigner- 
Weisskopf approximation, we find^ 

wr..^ r -Y- l(l|l^/(0)lx-p,^p)P ..;, ^r 

WoitA) = lim « > " ^ , , ^ , , = zAi?$ H — (111.9) 

^'^ .^0+ ^ m^ - E^{p) - E^{p) + ie 2 ^ ' 

where 

^E^^vj2 ° ' ;f ^ \ , (III.IO) 

^ m$ - E^{p) - E^[p) 

is the second order shift in the energy which will be absorbed into a renormalization of the 
$ mass and 

r ^ 27r^ |(l||i^,(0)|X-p,^p)r5("^* - E^{p) - EM) (III-ll) 

p 

is the decay width as per Fermi's Golden rule. Therefore in this approximation, we arrive 

at 

C$(t) = e-'^^**e-?*, (111.12) 

where we now consider a $ with k = (decay at rest) and find 



leading to 

r = 2nJ2\MM\'S{m^ - E^ - E^p)) = -^-^ (111.14) 



2 ^* 
2" 



^ The long time limit in the Wigner-Weisskopf approximation is equivalent to the Breit-Wigner approxi- 



mation of a resonant propagator [33| and holds for t ^ l/?7i$ 



where 



P 



2r7id 



4, 4, 4 o22 o22 o22 

m$ + m + m^ - 2ma,m^ - 2ma,m,,, - Zm^m^ 



6^,76^ 



^$"^V 



^x''^ 



1/2 



Inserting 0111.12^ into f lll.SP we find the quantum state 



I^W) 



-iAE^t 



(111.15) 



e-^*|l?;0^;0^) + 5^C(p;t)|x-p)|^p)|0$) (III.16) 



where 



C(p;t) = A^$(p) 



1 _ g-i(m<i,,fl-Sx(p)-'Ev,(p)-*r/2)t 



(^£;^(p) + E^{p) - m^^R + zr/2 



(III.17) 



where m$^/j = m$ + ^Eq> is the renormahzed mass of the $ particle. In what follows we 
drop the subscript R and always refer to the renormahzed mass. 

At this stage we can make contact with the momentum entanglement discussion in ref . |29| 
and at the same time exhibit the true non-perturbative nature of the results above by 
considering the state obtained using naive perturbation theory. 

Taking the initial state at t = to be |lf ), then to leading order in g, the time evolved 
state is given by 



\nt)) 



1 



t 

I I e'""'' Hie-'^'^'' dt' + ■■■ 
'o 



1?) 
n / 



(III.18) 



Introducing a resolution of the identity 1 = YIk I'^)('^I ^^ ^^"^ 



i^w) = |i|) + E 



MM 



E^ + E^ -m$ 



X-p)l^p)|0$) + 



(III.19) 



In this expression we have neglected a disconnected three particle intermediate state in 
which the initial |1?) remains and the interaction creates an intermediate state with three 
particles. This contribution is truly perturbative, does not contain resonant denominators 
as fllll.lQP and corresponds to the disconnected diagram in fig. ([2]). 

The probability of finding the daughter states is given by the familiar result (Fermi's 
Golden rule) 



vit) = J2\MM\ 



'sm{{m^-EM-E4p))t/2) 
{m^- E^ - E^) 



Tt 



(III.20) 



where F is given by (]III.14p . This result is obviously only valid for t ^ 1/F. It is now 
clear that the generalized Wigner-Weisskopf method that yields the state (1III.16P with the 
coefficients given by (1III.17P is truly non-perturbative. 

The momentum entanglement between the daughter particles is akin to that discussed 
perturbatively in ref.|29| with some important differences. In ref. [29| momentum entangle- 
ment was studied for the vacuum wave function. The corresponding contributions are truly 
perturbative and do not feature the resonant denominators that lead to secular growth in 
time and are similar to the contributions that we neglect and are described by fig. ([2]). 
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The quantum state (1111.160 describes an entangled state of the parent and daughter 
particles. The full (pure state) density matrix is given by 

p{t) = mt)){nt)\ (111.21) 

and its trace is given by 



TTp{t) = e~'' + V / ^|C(p;t)|^ (III.22) 



The momentum integral in (1111.22^ can be computed by changing variables to £ = E^{p) + 



E^ip). In the narrow width limit F ^ m<^,in^ + m^, the integrand is sharply peaked at 
£ = m$ so that the lower limit can be consistently taken to — oo thus allowing the integral 
can be computed by contour integration. We find 

^/(|^|C(p;t)r = l-e-^S (III.23) 

confirming that 

Trp{t) = 1 , (III.24) 

consistent with unitary time evolution and the unitarity relation (III. 151) . Furthermore the 
average number of ip (or x) particles is given by 

n^{p;t) = {ntH{p)a4p)\nt)) = (2^)'^^ = \C{p;t)\\ (III.25) 

thus the total number of tp (or x) particles is 

N^{t) = vj^s^^iP-^^) = 1 -e""*- ("1-26) 

Tracing out one of the daughter particles, for example x if it is unobservable, leads to a 
reduced density matrix 

p^t) = Tr,p(t) = e-r* |15)(1?| + J2 \Cip-,t)\' |^p)(^p| . (III.27) 

p 

This reduced density matrix describes a statistical mixture of states, the Von-Neumann 
entanglement entropy is therefore given by 

S{t) = -n*(0, t) ln[n*(0, t)] -^n^{p-t) In [n^{p ; t)] . (III.28) 

p 

where n^{p) is given by (IIII.25P and n*(0,t) = e~'"*. Because in the narrow width limit 
|C(p;t)p is a sharply peaked distribution, under integration with functions that vary 
smoothly near p ~ p* it can be replaced by 



Using this approximation we find 

S{t) = Tte-^' - [1 - e-^*] In [l - e"^*] - [l - e"^*] In [n^{p* ; oo)] (111.30) 

where 

This result has the following interpretation. In the asymptotic limit t ^ l/F, the entan- 
glement (von Neumann) entropy approaches (minus) the logarithm of the available states. 
The decay of the parent particle produces entangled pairs in which each member features a 
very narrow distribution centered at p* of width ~ F and height ~ l/F. The total area in 
momentum space yields 1/V since there is only one particle (of either type) produced in the 
volume V. Within the range of momenta centered at p* and of width F all of the available 
single particle states have equal probability oc l/V, therefore these states are maximally en- 
tangled as Bell states. This observation becomes clearer recognizing that a typical quantum 
state that contributes to the sum in flIII.16P is of the form 



c(p*,t)[ix(r))l^(-r)) + lx(-r))l^(r))J ; r = p*n, (111.32) 

where n is the direction of emission of either member of the pair. The quantum states 
with momenta p* — F/2 < p < p* + F/2 are represented with nearly the same probability 
|C(p*,t)p oc 1/yF in the sum. These states are Bell-type states and are maximally entan- 
gled, in fact these are similar, up to an overall normalization, to the entangled B^B^ states 
resulting from the decay of the T(45') resonance [ol-ll l| . but with the opposite relative sign 



because of charge conjugation. If the decaying particle has a short lifetime corresponding to 
a broad resonance, the emitted pairs will feature a distribution of momenta with probabilities 
|C(p, t)| determined by the Lorentzian profile of the resonance. 

As discussed above n'^(p* ; oo) (see eqn. ( 1III.25I) ) is the asymptotic phase space density 



of the produced particle (either x or ip). The entanglement entropy vanishes at the initial 
time since the density matrix at t = is a pure state and grows to its asymptotic value on 
a time scale l/F. 

The volume factor requires further discussion. It appears as a result of treating the 
initial (and final) states as plane waves. A more physical description of the decay process 
would treat the initial state as a wave packet localized in space. A natural localization 
length scale for the decaying state would be the decay length r ~ l/F; for example, for the 
experimentally relevant case of a pion decaying at rest, r c^ 7.8 mts is also approximately 
the typical experimental scale of the production region. In particular, for a pion decaying 
into a muon and a neutrino with p* = SOMeV and m,y ^ m^ and V ~ 1/F^ ,we find 
n^'i^p* ; oo) ~ 10 



-30 



IV. CONCLUSIONS AND FURTHER QUESTIONS 

The decay of a parent particle leads to a quantum entangled state of the daughter particles 
as a consequence of conservations laws in the decay process. We generahzed and extended a 
method used in the study of spontaneous decay of atomic systems to the realm of quantum 
field theory to obtain in a consistent approximation, the full quantum state that describes 
the time evolution of the decaying particle and the production of the daughter particles. This 

10 



method is non-perturbative and is manifestly unitary. We have implemented the method 
to study the simpler case of bosonic parent and daughter particles to highlight the main 
concepts and consequences, however, the results are quite general. 

The full quantum state resulting from the time evolution of the decaying particle yields 
a pure state density matrix. However, if one (or more) daughter particles is unobserved, 
tracing over their degrees of freedom leads to a mixed state density matrix whose time 
evolution is completely determined by the unitary time evolution of the decay process. This 
mixed state density matrix features an entanglement entropy which is a manifestation of 
the quantum correlations of the entangled product state. We obtained the time evolution 
of the entanglement entropy and show that it grows on a time scale determined by the 
lifetime of the decaying particle and reaches a maximum that corresponds to the logarithm 
of the available phase space states of the decay particles. For a parent particle described by 
a narrow resonance the distribution of produced (entangled) daughter particles is nearly a 
constant in a narrow energy-momentum region, the emitted particles are nearly maximally 
entangled Bell-states. These concepts may prove relevant in the statistical analysis of the 
time evolution of entangled B-meson pairs in studies of CP violation. 

This work can be viewed as prelude to the calculation of the entanglement entropy due 
to particle decays in de Sitter space [331]. There, due to the fact that particle can decay into 
itself with momenta that are much less than the Hubble constant -ffdeSitter [33] , we expect 
there to be an interplay between the horizon size and the decay rate that will feed in to the 
behavior of the entanglement entropy. This may then mix the ideas of entanglement entropy 
developed here with those coming from studies of spatially separated portions of de Sitter 



space such as in ref.J28|. Work on these aspects will be reported elsewhere [35 
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